Space-time resolved electrokinetics in cylindrical and semi-cylindrical microchannels 



Michele Campisi 
Department of Physics, University of North Texas, 
P.O. Box 311427, Denton, TX 76203-U27,USA 
(Dated: February 1, 2008) 

It is shown show how to employ Bessel-Fourier series in order to obtain a complete space-time 
resolved description of electrokinetic phenomena in cylindrical and semi-cylindrical microfluidic 
channels. 
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PACS numbers: 47.15.-x, 47.65.+a, 73.30.+y, 82.39.Wj. 



The employment of Fourier series and Laplace trans- 
forms allows to solve exactly the low Reynolds num- 
ber Navier-Stokes equation in rectangular microchan- 
nels [H, In Refs. [H, 01 this has been used to find 
the space-time resolved velocity field and current den- 
sity generated by the sudden application of a possi- 
bly time dependent pressure gradient and/or an elec- 
tric field. This allowed to fully characterize transient 
and steady state electrokinetic phenomena in rectan- 
gular microchannels. A similar characterization in the 
case of cylindrical and semi-cylindrical is still missing. 
These channel geometries (particularly the latter) are 
relevant for microfluidics devices fabricated from glass 
substrates. Typically microchannels are obtained from 
such substrates by means of isotropic etching which pro- 
duces semi-cylindrical grooves. Here we shall show how 
to adapt the mathematical methods presented in Ref. 
P, [|| to the cylindrical and semi-cylindrical geometries, 
in order to allow for a space-time resolved description 
of electrokinetic phenomena in these cases. The physics 
involved will be qualitatively the same as in the rectan- 
gular geometry, thus we will not delve into a detailed 
discussion, but present only the mathematical method. 

Cylindrical geometry As reported before [l| the 
starting point for studying space-time resolved electroki- 
netic phenomena in microchannels is the low Reynolds 
number incompressible Navier-Stokes equation with a 
generic time and space dependent body force. The nat- 
ural set of coordinates for solving this problem within 
cylindrical microchannels is the cylindrical coordinates 
(r, <p, x), where x runs along the channel axis and r,tp 
are polar coordinates in the channel cross-section. Using 
this coordinate set, the x component of the Navier-Stokes 
equation takes the following form 



solution of the problem with null boundary condition over 
a rectangular boundary. Likewise, for a circular bound- 
ary of radius R the appropriate basis set is given by the 
following functions 



ip m ,k(r, if) 



'Jm{oim,kr/R) m e Z, k e N+ 



(2) 



The symbol J m denotes the m th Bessel function of the 
first kind, and a m ^ denotes its k th zero. For sake of 
completeness let us recall that the J m are defined as the 
solutions of the equations 

p2 <V Jm ^ + 9 dp Jm ^ + ^ 2 ~ ^ ^ Jm ^ = ° ^ 

and that for fixed m the following orthogonality relation 
exists between the functions J m : 



dppJ m (a m ,kp)J m (a m , q p) = -5k, q [J m {um,k)] 2 (4) 



In a similar fashion the y m {p>) — e imv are solutions of 
y' = imy and obey the orthogonality condition: 



d<Py™(v)yn(¥>) = 27175,, 



(•5) 



where the symbol "*" denotes complex conjugation. Us- 
ing Eq. (|3]) and Eq. ((SJ), allows to normalize the basis 
set © and obtain the the following complete orthonor- 
mal basis: 

^,fc(r,y) = e Tp7 (aw ;* r/ ? meZ,fceN + (6) 

y/TTHJ m +l{a m ,k) 



du{r,ip,t) 
dt 



3P_ 

dr 2 



ld_ 

r dr 



r 2 dtp 2 



u(r,ip,t) = f(r,<p,t) 
(1) 



where u is the x component of the velocity field and / 
is the body force. In the case of rectangular geometry 
the solution of the Navier-Stokes equation was found by 
means of double Fourier series expansion. In other words 

the set of functions ip^q = sin O^f) sin (ff^), k, q € N+ 
was chosen as the appropriate complete basis set for the 
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dip rdri>* k {r,ip)i) niq (r,tp) = 5 mtn 5 kA (7) 



The completeness of the set ([5]) allows to expand the 
functions / and u in a double generalized Fourier series 
(the Bessel-Fourier series), as follows: 



u(r, ip, t) 



EE 



(t)ip m ,k(r, <p) 



(8) 



m 



with coefficients given by: 

p2ir pR 

u m . k {t) = dip rdrip* k (r,(p)u(r,ip,t) (9) 



Expanding the Navier-Stokes equation over the ba- 
sis © gives, thanks to the property ([3]) and in complete 
analogy with what was found previously for the rectan- 
gular geometry the following set of equations: 

d 

7-"miW + ^ A, 2 „ k U m . k (t) = f m , k {t) (10) 

at 

with null-boundary condition over the r = R circumfer- 
ence automatically fulfilled. The quantities A^ fe , given 
by the following formula, 



A h = 



R 2 



(11) 



are, so to speak, the expansion coefficients of the Lapla- 
cian operator over the given basis set ([6]). The f m ,k(t) 
represent the expansion coefficients of the body force /: 

p2ir pR 

f m ,k{ t )= I dl P I rdrip^ k (r,(p)f{r,ip,t) (12) 
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If the liquid is considered as initially at rest Eq. (fTOf must 
be solved with for the initial conditions u mjk (t = 0) = 0. 
Following the general line drawn in Ref. [l[ , the solution 
is easily expressed in the Laplace space as: 



k(s) 



fm,k( s ) 



(13) 



which expresses the Laplace transform of each coef- 
ficient of the velocity profile in terms of the corre- 
sponding Laplace transformed components of the driving 
force. The solution u(r, ip, t) is obtained by anti-Laplace- 
transform (C^ 1 ) and summing up: 



u(r, ip, t) 



keN + 



fm,k{ S ) 



i>m,k{r,<p) (14) 



For a pressure driven flow the body force / would be 
given by f(r,p,t) = , with L the length of the 

channel, p the liquid density and AP(t) the possibly time 
dependent pressure difference applied at the ends of the 
channel. For an electro-osmotically driven flow the body 
force would be given by p '^'^ 5(t), with E(t) the ap- 
plied electric field and p e the electric double layer (EDL) 
charge density that spontaneously forms at the solid liq- 
uid interface. The latter can be found by solving the 
Poisson-Boltzmann equation [1] for the electric double 
layer potential tp within the Debye-Hiickel approxima- 
tion [a]: 



where x is the inverse Debye length. Expanding over the 
basis (El, and using the Poisson equation Aip = —p e /e, 
like in [lj one obtains the cherge density coefficients: 



Pe(m,k) 



A m ,fc-fm,fc 

A 2 U + Y 2 
m.k ' A 



(16) 



where I mm k denote the expansion coefficients of the unity: 

P 27T n Ft 

I m ,k(t)= dip rdrip* mk {r,tp). (17) 
Jo Jo 

The solution of the problem in the cylindrical geometry is 
formally equivalent to that of the rectangular geometry. 
The only difference is contained in the way the expansion 
cocfficeint are calculated: using the double Fourier series 
for the rectangular case, and using the Bessel-Fourier se- 
ries in the cylindrical case. Thus once the basis set ap- 
propriate for a given geometry is found, the problem is 
automatically solved. 

Semi-cylindrical geometry In this case, the func- 
tion u(r, if) must obey not only the condition of being 
null for r = R, but also for ip = 0, n. Seen in a differ- 
ent way, the function u must be odd under the operation 
ip — > — ip. Therefore its expansion series would contain 
only odd terms, i.e., it would be of the type: 

X] [^m,k{r,ip) -ip- m ,k{r,ip)] (18) 

Namely it would contain only sine terms. Therefore the 
orthonormal basis set suitable for the semi-cylindrical ge- 
ometry is: 

2sm{rrup)J m (a m , k r/R) 

<?W(r, ip) = - F r— r m, k G N+ (19) 

^/nHJ m+ i(a m ^k) 

Where the cj) m k satisfy the orthonormality conditions: 



dip I rdr(f> m}k (r,(p)(f> n)q (r,ip) = S mjn 5 k:q (20) 
o Jo 

We will write the expansion of u as: 

u(r,ip,t)= ^ u m,fe(*)^m,fc(r,v) (21) 

with coefficients given by: 

/>7T pR 

u 'm 1 k{ t )=j d P rdr(j> m ,k(r,<p)u(r,(p,t) (22) 
Jo Jo 

where the prime symbol is used to distinguish these co- 
efficients from those defined previously. Adopting the 
same notation for the expansion of the body force, again 
the expansion of the Navier-Stokes equation leads to the 
solution 



3 



which is formally equivalent to Eq. 
density will be given by 



T3}. The charge 



Pe(m,k) 



A 2 J' u 

2 m,/c m,k 



(24) 
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where 



dip / rdrcf> m< k{r,ip) 



(25) 



which is formally equivalent to Eq. (|16|) . As an illus- 
tration of the method Fig. [1] shows a typical plot of 
the EDL charge density obtained from Eq. All 
the information relevant for the description of electroki- 
netic phenomena in cylindrical and semi-cylindrical mi- 
crochanels is contained in the coefficients u m< k, p e (m,k) 
and u' j.,pff m M , respectively. These can be used like 
in Ref. [l|, [2J to obtain a space-time resolved description 
of electrokinetic phenomena. For example, for the semi- 
cylindrical geometry one finds the following generalized 
conductance matrix: 



m = — y 

m,k 



-A' 



(26) 



Where L is the channels's length, uo is the angular fre- 
quency of the driving body force, 
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and 



A' 



T/2 

m,k 



IL 



L m,kPe{k,q) 

T n' n n 

1 m,kl J e(k,q) > J e(k,q) 



(27) 



FIG. 1: Typical EDL charge density profile in non- 
dimensional units (R = l,p e (l,0) = 1). The first 100 x 100 
Bessel-Fourier coefficients have been employed to generate the 
plot. 



CfcH = arctan 



(28) 



